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Abstract 

We investigate the discrepancy that exists at 1ow-xt = 2pt/\/s be- 
tween the next-to-leading order QCD calculations of prompt photon pro- 
duction and the measured cross section. The central values of the mea- 
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■ sured cross section are of order 100% larger than QCD predictions in this 

region. It has been suggested that the bremsstrahlung contribution may 
^ ' account for this discrepancy. The quark fragmentation function Dj/g{z) 

' has not been measured and an exactly known asymptotic form is normally 

used in calculations. We examine the effect of much larger fragmenta- 
tion functions on the QCD predictions. After illustrating the effect of 
the large fragmentation functions in some detail for recent CDF data at 
^^5=1. 8 TeV, we perform a fit to 8 prompt photon data sets ranging 
in CMS energy from 24 GeV to 1.8 TeV. While a large fragmentation 
function normalization may prove to play an important role in resolving 
the discrepancy, the present theoretical and experimental uncertainties 
prevent any definite normalization value from being determined. 
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1 Introduction 



Prompt photon production in hadronic collisions is an important testing ground 
for QCD. Recently, measurements of inclusive prompt photon production at 
CMS energy ^/s =1.8 TeV have been published by the CDF collaboration, and 
a comparison with QCD calculations at next-to-leading order has been carried 
out llj]. While the QCD predictions give good qualitative agreement, the central 
values of the measured cross section in the low-pr region are of order 100% 
larger than standard QCD calculations. The same discrepancy exists between 
theory and experiment for data taken by UAl and UA2 at the CMS energy 
v/i=630 GeV. 

The NLO QCD calculation consists of three parts. 

^ ^LL I ^NLL I ANOM 

where the three terms indicated are the leading logarithm (improved Born ap- 
proximation) contribution, the next-to-leading logarithm piece, and the so-called 
"anomalous" part, in which a photon is emitted coUinearly from an outgoing 
parton. This is just the bremsstrahlung process. While the first two parts 
of the QCD calculation are reliably calculated (given the fact that the quark 
and gluon distribution functions are reasonably well constrained) the third, the 
bremsstrahlung process, is not. The calculation of the bremsstrahlung process 
involves an unmeasured set of non-perturbative functions, the parton — > photon 
fragmentation functions. While there is a standard parametrization for these 
functions used by most theorists in QCD calculations, this parametrization is 
only valid at asymptotically large energy scales where the effects of boundary 
conditions are vanishing. In this paper we examine the effect of large fragmen- 
tation functions on the QCD predictions of inclusive prompt photon production 
and compare with experimental data. 

This paper is organized as follows. In Section 2 we discuss the NLO QCD 
calculation of the bremsstrahlung component of prompt photon production and 
we elaborate on certain aspects of the fragmentation functions involved in this 
calculation. In Section 3 we first illustrate the existing discrepancies between 
theory and experiment for the UAl and UA2 data at A/i=630 GeV and CDF 
data at 1/5=1.8 TeV. Concentrating on the CDF data, we examine the effects 
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of the theoretical uncertainties due to the choice of scales and distribution func- 
tions on the discrepancy. We then go on to show the results of the QCD prompt 
photon calculation including larger fragmentation functions. We compare our 
results with the standard calculations, and with the CDF data. We again ex- 
amine the effect of changing the scales and parton densities on the theoretical 
predictions. In Section 4 we proceed to compare the results of the calculation 
with and without a large set of fragmentation functions to data at center of 
mass energies yi=1800, 630, 546, 63, and 24 GeV. A analysis is performed 
and the results are discussed. Section 5 is reserved for the conclusions. 

2 Bremsstrahlung calculation; Fragmentation functions 

The anomalous part of the QCD calculation of the inclusive prompt photon 
production differential cross section a{AB ^ 7 + X) is given at the leading 
logarithm level by 

^ANOM = J2 J Fa/AiXa, Fb/siXb, Dj/c{z, Mf) ^ dXa dXb dz dv. (1) 
abed 

In this equation v = 1 +t/s, where t and s are the Mandlestam variables of the 
subprocess, and Xa, x^, and z are the fractional momenta of the initial state par- 
tons and the photon, respectively. The subprocess cross section da/dv, which 
describes the parton level scattering process ab cd, is convoluted with the 
parton distribution functions Fa/ a and F^/b and the parton photon fragmen- 
tation functions D^/^- The three scales in the calculation are the renormalization 
scale fiR, associated with the strong coupling constant, the factorization scale 
fip where the parton distribution functions are evaluated, and the fragmenta- 
tion scale, Mf, associated with the parton — > photon fragmentation process. At 
sufficiently high order in perturbation theory, the total calculated cross section 
should be insensitive to the chosen values of these scales. 

The fragmentation and distribution functions are nonperturbative objects 
and are treated as measured inputs within the framework of perturbative QCDQ 
While many experiments such as deep inelastic scattering have given rise to a 

* Given a set of functions at some scale the evolution of them to other energy scales is 
perturbatively calculable. 
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fairly well constrained set of quark and gluon distribution functions, the par- 
ton photon fragmentation functions have not been measured. However, ow- 
ing to the pointlike nature of the photon these fragmentation functions differ 
essentially from the parton distribution functions in that at asymptotically large 
scales the fragmentation functions are perturbatively calculable. The asymptotic 
form for the quark and gluon fragmentation functions, which includes the sum- 
mation to all orders in perturbation theory of soft and coUinear gluon emission 
at leading logarithm level, is parametrized as follows [0] 



e2(2.21 - 1.282 + 1.292 
1 - 1.63 ln(l - z) 



2n^0.049 



+ 0.0020(1 - zy-^z 



2.0^-1.54 



(2) 



.0.194 



zD!;;^^{z,Q')=F^-^{l-zrh-'-'\ (3) 

where F = (a/27r) In (^Q^/Ag^^^. Theorists generally use this asymptotic form 
when doing QCD calculations of prompt photon production and comparing to 
the experimental data p|, |^, However, the boundary conditions on the evo- 
lution equations may be important at scales where experimental data has been 
taken. The boundary conditions on the evolution equations are simply the ini- 
tial fragmentation functions at some scale Qq, where Qq is defined to be the 
scale above which the perturbative evolution equations are assumed to be valid. 
We take this scale to be 2 GeV. We consider the possibility that the physics at 
scales below Qo builds up a set of functions at Qo with a larger normalization 
than that of the asymptotic form evaluated at the starting scale. 

In order to make the analysis tractable we consider starting quark and gluon 
fragmentation functions which have the same shape as the asymptotic forms 
but which have normalizations A^o times that of the asymptotic forms evaluated 
at the scale Qq. Of course, as the scale increases these large fragmentation 
functions evolve and approach the asymptotic form normalization. In figure la 
we show the up-quark fragmentation function both for the asymptotic form and 
for a large fragmentation function with Nq = 20. The asymptotic form grows 
logarithmically with the scale and this factor has been divided out in fig. la. 
We show in this plot the large fragmentation function evaluated at the starting 
scale 2 GeV, the scale 10 GeV and the scale 100 GeV. In figure lb we show the 
normalization of the large fragmentation function vs. the scale Q, for values of 
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Nq equal to 10, 20, and 30. Because of isolation cuts, many experiments are 
only sensitive to large values of fragmentation function fractional momenta. For 
this reason, the relative normalization N{Q'^) shown in fig. lb is given by 



where the denominator is the asymptotic form fragmentation function of eq.2. 
Note that the relative normalization decreases rapidly in the first 10 or 20 GeV 
beyond the starting scale Qq. This feature is well suited for resolving the dis- 
crepancy between the QCD prompt photon predictions and experimental data 
at low transverse momenta. If the starting fragmentation function normalization 
is much larger than the asymptotic form normalization, the smaller px values 
will probe this large normalization, while at larger px values most of the initial 
large normalization will have evolved away. Thus, in the low-p-r region the pre- 
dicted cross section will markedly increase relative to the standard calculations, 
while at higher photon transverse momenta we can expect results similar to the 
standard QCD predictions. 

It is of interest to estimate the largest initial fragmentation function nor- 
malization that may be physically relevant. An order of magnitude estimate 
for an upper bound comes from demanding that at the starting scale Qq the 
total number of photons fragmented from some parent quark with fractional 
momenta above some cut-off zq be less than some small number of 0{a/as)-, say 
^, times the total number of hadrons that are fragmented by the quark. We 
show in figure 2 the maximum normalization N^ax as a function of Zq coming 
from this constraint. We assume a hadron multiplicity of \/2E for a parent 
quark of CMS energy E, in GeV units, and have fixed the parent quark energy 
at 3 GeV. Taking zq to be 0.3, we see from fig. 2 that we can tolerate an initial 
normalization which is ~20 times larger than that of the asymptotic form. 

Another more quantitative bound can be estimated by examining the par- 
ton photon process in the crossed channel, namely the photon quark 
fragmentation function. In this case the photon structure function in the non- 
perturbative regime can estimated by appealing to vector meson dominance 
(VMD) arguments. In ref. [^] the 7 — > vector meson fragmentation function is 
hypothesized to be equal to the 7 — > 7r° fragmentation function, which is ex- 
tracted from fits to e^e~ vr^X data. We find from the analysis of ref.^ that 
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using the VMD hypothesis to specify the fragmentation function boundary con- 
ditions hmits A^o to be less than 2 at leading order. However, due to the many 
theoretical uncertainties involved in the VMD procedure, we will not respect 
this limit. Rather, we will go ahead and vary Nq freely in order to determine 
whether the present experimental prompt photon data can constrain it directly. 

We will find that in order to account for a substantial increase in the prompt 
photon production cross section at low transverse momenta we will need to con- 
sider initial fragmentation function normalizations which are in some cases 30 or 
more times larger than the normalization of the asymptotic form (which is proba- 
bly unreasonably large). At first sight this may be surprising since at large center 
of mass energies and low-pr it is well known that the bremsstrahlung contribu- 
tion can dominate the other photon production processes. However, in order to 
reduce the background photons from 7r° and f] decay the data samples are taken 
with the imposition of an isolation cut. These isolation cuts effectively remove 
most of the bremsstrahlung events from the data sample. Typically only photon 
fractional momenta values of z>0.85 are allowed. Additionally, at the lowest pt 
values measured in the experiments the fragmentation function normalizations 
will have evolved considerably down from any large starting normalization. 

3 Results 

We use the algorithm of Baer et al. Q to calculate the next-to-leading order 
QCD cross section for prompt photon production in hadronic collisions. See 
their paper and references therein for the relevant formulae. The calculation 
involves a Monte Carlo integration wherein various isolation criteria are easily 
implemented. To simplify our analysis we choose to set the renormalization 
scale equal to the factorization scale, ^'^ = = 

In figure 3 we compare our results for prompt photon production with 
CDF data [0] taken at ^5=1.8 TeV, and with UAl and UA2 data 0] taken at 
630 GeV. Here we are using the standard asymptotic form parametrization of the 
quark and gluon fragmentation functions given in eqs.(2) and (3) {i.e. Nq=1). 
We show the results for the scales /x^ = M| = and the distribution functions 
Morfin & Tung Set 1 (M-Tl) P). Both the CDF data and the corresponding 
calculation include an isolation cut, such that the total hadronic energy inside a 
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cone of size AR = ^/At]'^ + A0^=O.7 centered on the photon be less than 2 GeV 
{■)] and (p are the psuedo-rapidity and azimuthal angle, respectively). The UAl 
data and the corresponding calculation include an isolation cut such that the 
total hadronic energy inside a cone around the photon of size AR=OA be less 
than 10% of the photonic energy. The isolation criteria for the UA2 data shown 
in fig. 3 is not as simple, but in an earlier UA2 publication of prompt photon 
data 1^ it was demonstrated that the data sample was consistent with having 
assumed an isolation cut such that the hadronic energy in a cone around the 
photon of size Ai?=0.78 be less than 3.5 GeV. We use this isolation criteria in 
our calculation. The UAl and UA2 isolation criteria give very similar results 
and we show the curve corresponding to the UAl isolation criteria in fig. 3. 
The data samples have the overall scale normalization uncertainties shown in 
the figures. The error bars in the figures include statistical and systematic errors 
added in quadrature. The discrepancy in the low-px region is apparent for all 
three data sets. 

The discrepancy is more easily seen on a linear scale. In figure 4a we show 
the relative difference between data and theory, (data-theory) /theory, versus 
Pt for the CDF data. This figure illustrates the theoretical uncertainty due to 
changing the scales and distribution functions. The reference theory calculation 
on this plot includes M-Tl distribution functions and scales equal to p"^. The 
results using the DFLM260 [|10| distribution functions are shown and give about 
a 20% larger cross section at low-px, decreasing to less than 10% of the M- 
Tl result at large pt- Also in this figure we see that changing the scales to 
/i^ = Mj = p'^/4 increases the results of the predictions by 10-15% (for a given 
set of distribution functions). The results are relatively insensitive to changes 
in the factorization scale Mj. For standard initial normalizations {Nq = 1) 
the results change by < 5% when considering Mj = p^ and Mj = or 
Mj = 4* p\,. At larger initial fragmentation function normalization values the 
fragmentation scale dependence is larger, but generally it is < 10% for starting 
normalization values less than 40. 

The discrepancy is readily observed in fig. 4a for the case of scales 
and M-Tl distribution functions. For these choices the central values of the 
three lowest pt bins are 60, 51 and 117% larger than the theoretical calculation, 
and have a x^/d.o.i. of 5.3. For comparison, the central values of the next 



6 



lowest three pt bins average less than 9% larger than the theoretical curve. The 
choices corresponding to the largest theoretical curve arep|./4 scales and DFLM 
distribution functions. In this case there is certainly no glaring discrepancy in 
the low-pT region, although the central value of the lowest pt bin is still 80% 
larger than the prediction. 

In figure 4b we illustrate the discrepancy between data and theory for the 
UAl and UA2 data at -^/i^GSO GeV. As in fig.4a the largest low-pr region dis- 
crepancy corresponds to the theoretical calculation which includes M-Tl parton 
densities and p^ scales. For these choices the lowest four pt bins average 97% 
larger than the expectation. The discrepancy is least, as for the CDF data, in 
the case of DFLM parton densities and scales equal to p\/^. In this case the 
X^/d.o.f. of the lowest four pt bins is 3.2. 

We show in figure 5a the qualitative improvement in the comparison of 
CDF data and theory if we include a large initial normalization for the fragmen- 
tation functions. We show the curves for A^o=30 and 40 for the case of M-Tl 
distribution functions and p\ scales. The large A^o curves qualitatively account 
for the rise in the cross section in the low-pT region while retaining the good 
description at large px- The best fit (defined by the smallest x^; see Sec.4) for 
the initial normalization Nq for this choice of scales and distribution functions 
is Ao=34. Of course, the best fit value for the starting normalization Nq will 
change with the distribution function and scale choices. If the scale /i^ is chosen 
to be smaller [e.g. [x^ = Pr/'^) the QCD prediction will increase overall, thus 
a smaller normalization will be favored. As an illustration, we show in figure 
5b the standard (A'o=l) curve and the best fit {Nq = 17) curve for the scale 
choice = Mj = p'^/4:. Similarly, as the DFLM distribution functions give 
larger cross sections, they will also give smaller values of the best fit starting 
normalization. Figure 5c shows the DFLM curve for the standard calculation 
[Nq = 1) and the best fit initial normalization in this case, A^o = 17, for the 
choice of scales /i^ = Mj = p\. 

If we move from the choice of distribution function and scale of fig. 5a to 
the one of DFLM distribution functions and the scale choice /x^ = Py/4, we 
can expect an even smaller best fit initial fragmentation function normalization 
value, as both of these choices increase the theoretical values of the cross section 
relative to the choices of fig. 5a. We find the best fit value in this case is A^o = 4. 
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However, it is perhaps worth noting that this case gives the worst fit to the data. 
The central values of all the data points (except for the three lowest pt bins) fall 
below the theoretical calculation (even in the case A^o — 1) and correspondingly 
the value for this choice of scales and distribution functions is the largest 
of the four choices considered. 

We have seen that while varying the choice of scales and distribution func- 
tions in the theoretical calculation at TeV the best fit value for the initial 
normalization A^o varies from 4 to 34. The uncertainty in the QCD calculation is 
too great to yield a well constrained value for the normalization. However, even 
if the QCD calculation were under better control, we will see in the next section 
that the experimental errors are too large to lend a best fit normalization value 
any statistical significance. 



4 analysis 

In order to understand the statistical significance of the qualitative improvement 
achieved by stipulating a large initial fragmentation function normalization, and 
to see if such a large fragmentation function is consistent with other lower energy 
experimental data, we vary the initial fragmentation function normalization and 
perform a x^ analysis. In most of the data sets there is an overall scale nor- 
malization uncertainty. In such cases we evaluate the x^/d.o.f. by varying the 
overall normalization of all the data points and finding the normalization A'^ 
which renders the x^ formula 

a mininnim. AN is the experimental normalization error. For Acr^^^ we add 
the statistical and systematic errors in quadrature. If part of the systematic 
error is correlated the true x^ values for that data sample will be larger than 
the ones evaluated using eq.(5). In figures 6 we show the distributions for 
three representative data sets; the CDF data at ^/s =1.8 TcV, the combined x^ 
for the UAl and UA2 data sets at =630 GeV, and the NA24 data taken 
at ^/s =23.76 GeV. In fig. 6a the scales are set to /i^ = M| = p\ and the 
DFLM distribution functions are used. These three data sets show a favored 
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initial fragmentation function normalization value much larger than the standard 
value of 1 (the values are A^o=17, 24, and 19 respectively). However, for both the 
1.8 TeV and the 630 GeV data, the difference in the value at the minimum 
and that for an initial normalization of 1 (or zero) is less than 0.5. Thus, while 
the experimental data favors a large initial fragmentation normalization, the 
experimental errors are such that a standard initial normalization (or even zero 
initial normalization) gives almost as good a fit. 

For the low energy data set shown, the situation is quite different. The 
reason the distribution corresponding to the NA24 data set shows a strong 
dependence on the initial fragmentation function normalization is that this data 
set has no isolation criteria applied. The plot in fig. 6a shows the value A^o = 1 
is actually ruled out at the 68% confidence level. However, by changing the scales 
and/or parton distribution functions in the calculation the value Nq=1 can be 
accommodated. This is shown in fig.6b, in which the scales in the calculation 
have been set to pf-/4. All three distributions in fig.6b show minima at smaller 
values of A"o than in fig. 6a, and in particular the NA24 data in this case shows 
a minimum at A'o=5. The x^/d.o.f. value corresponding to A^o = 1 differs by 
less than 1 from the x^/d.o.f. value at the minimum. 

We now discuss the combined plots for all of the data samples considered. 
The data samples included are the CDF data at ^/s =1.8 TeV and = [Q; 
the UAl data at y/s =630 GeV and y/s =546 GeV, both for r/ = and 1.1 
0; the UA2 data at ^ =630 GeV, = §; R806 data at =63 GeV 
and = 0; and NA24 data at y/s =23.76 GeV with -.65 < r] < .52 



Eight QCD calculations were performed for each of the data sets corresponding 
to the distribution function choices M-Tl and DFLM, and the scale choices 
/i^ =n*pip] Mj = m*p'ip for (n, m) = (l,l), (1,4), {\, \), and (|, 1). We show in 
figure 7a the distribution as a function of the initial fragmentation function 
normalization Aq for both sets of structure functions considered and for the 
scale choice /i^ = M| = p^. In this figure the initial normalization values 
corresponding to the lowest x^ values are Ao=24 and Ao=19 for the M-Tl and 
DFLM distribution function choices, respectively. In figure 7b we show the 
results for the scales equal to Pr/4. Here the minimum x^ values correspond to 
the normalizations Aq = 8 and Aq = 4 for the M-Tl and DFLM distribution 
function choices, respectively. The distributions corresponding to the same 
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//^ scale choices as in figs. 7 but different M| scale choices give results similar to 
those shown in figs. 7. 

The plots all show very mild dependence on the initial normahzation A^o- 
A mild dependence for the higher energy data sets can be explained as follows. 
In the case of the CDF data the theoretical predictions and the experimental 
data agree within the (large) experimental errors for 11 out of 12 of the highest 
Pt data points. The data set, even with a discrepancy for the two or three lowest 
Pt data points, has a low x^/d.o.f. The x^/d.o.f. value changes very little if the 
theory is modified in such a way that the three lowest pr data points are also 
brought into agreement. 

In none of the eight cases considered is the difference between the x^ values 
corresponding to the standard normalization A^o =1 a-nd the normalization at 
the minimum of the x^ plot greater than 0.4. Thus, because the experimental 
errors are large, the quantitative improvement in the comparison of theory and 
data obtained by including a larger fragmentation function normalization is 
small. However, we do find it somewhat suggestive that in the individual x^ 
distributions for all of the data sets considered and for all choices of scales and 
parton densities considered (64 plots) the minima were more often than not at 
initial fragmentation normalization values greater than 16. 

5 Conclusions 

We find that including large boundary conditions in the evolution equations of 
the quark — > photon and gluon — > photon fragmentation functions may play an 
important role in resolving the discrepancy between next-to-leading order QCD 
calculations and experimental data of the inclusive prompt photon production 
cross section in the \ow-xt region. The fragmentation function evolutions were 
performed at leading order here. Perhaps including large initial fragmentation 
function normalizations in a full next-to-leading analysis would lead to a better 
constrained set of best fit initial normalization values. The inclusive prompt 
photon data sets considered here favor a larger initial normalization, but the 
theoretical uncertainties arising from scale variations and parton density un- 
certainties arc too large to give a definite value. A analysis reveals that, 
additionally, the experimental errors are too large to give any statistical sig- 
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nificance to what appears to be a qualitative improvement in the comparison 
of theory and data by including a larger fragmentation function normalization. 
In order to satisfactorily resolve the discrepancy, both theoreticians and experi- 
mentalists will have to further refine their respective results. On the theoretical 
side, more tightly constrained proton structure functions and higher order cal- 
culations would certainly stabilize the predictions. On the experimental side 
we note that if A^o is large the photon production cross section is particularly 
sensitive to the isolation criteria (either the isolation cone size or the amount of 
hadronic energy allowed in the cone) . A constraint on the initial fragmentation 
normalization could come from an experimental determination of the depen- 
dence of the photon production cross section on the isolation criteria. 
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Figure captions 
Figure 1 

(a) The up quark fragmentation function versus z. The curves are shown 
for A^o=20 at scales — 1 GeV, and Q—V^ and 100 GeV. The asymptotic form 
(A^o=l) is also shown. 

(b) The normalization of the up quark fragmentation function versus the 
scale Q for initial normahzation values A^o =10, 20 and 30. The normalization 
is A^o times the asymptotic form normahzation at the scale Qq — 2 GeV. The 
dotted lines show the asymptotes at GeV and A^=l. 

Figure 2 

The maximum initial normalization for the up quark fragmentation function 
Njnax versus the photon fractional momenta cut-off Zq, assuming the number of 
photons fragmented from a 3 GeV quark with fractional momenta z > Zq\& ^ 
times the number of hadrons fragmented by the quark. 

Figure 3 

The inclusive prompt photon differential cross section Ed^a/dp^ versus the pho- 
ton transverse momenta. The data is from CDF at TeV and UAl and 
UA2 at \/s =630 GeV. The curves are the QCD predictions using Morfin-Tung 
Set 1 distribution functions and scales equal to Note the experimental 
normalization uncertainties. 

Figure 4 

The relative difference between data and theory, (data-theory) /theory, is shown 
vs. transverse momenta. The reference theory calculation includes scales 
with M-Tl parton densities. Also shown are the QCD results for scales p^/A 
and DFLM distribution functions. 

(a) CDF data at ^5=1.8 TeV. 

(b) UAl and UA2 data, v^=630 GeV. 
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Figure 5 

The relative difference between data and theory, (data-theory) /theory, vs. trans- 
verse momenta for the CDF data at ^/s—1.8 TeV. The reference theory calcu- 
lations include the standard asymptotic form fragmentation functions. 

(a) M-Tl distribution functions and scales /i^ = Mj = are used. Initial 
fragmentation function normalizations A^o=30 and A'o=40 are shown. The best 
fit initial normalization value in this case is A^"o=34. 

(b) The results for M-Tl distribution functions with scales /i^ = Mj = p^/A. 
The best fit initial fragmentation function normalization value NQ=n is shown. 

(c) DFLM parton densities are used, and the scales are set equal to p\. The 
best fit normalization value iVo=17 is shown. 

Figure 6 

The distributions vs. the initial fragmentation function normalization Nq for 
three data sets; the CDF data at ^/s = 1.8 TeV, the combined for the UAl 
and UA2 data at ^/s =630 GeV, and NA24 data at = 23.76 GeV. All the 
data sets are at zero rapidity. The DFLM distribution functions were used. 

(a) The scales arc set to = M| = p\. 

(b) The scales are set to jj? = Mj =p'^/A. 

Figure 7 

The combined distributions vs. the initial fragmentation function normaliza- 
tion Nq for eight data sets; see text. 

(a) The scales arc set to ji^ = Mj = p^. 

(b) The scales are set to /^^ = Mj =p'^/A. 
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